Angular-dependent oscillations of the magnetoresistance in Bi2Se3 due to the 

three-dimensional bulk Fermi surface 
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We observed pronounced angular-dependent magnetoresistance (MR) oscillations in a high-quality 
Bi2Se3 single crystal with the carrier density of 5x10^* cm~'', which is a topological insulator 
with residual bulk carriers. We show that the observed angular-dependent oscillations can be well 
simulated by using the parameters obtained from the Shubnikov-de Haas oscillations, which clarifies 
that the oscillations are solely due to the bulk Fermi surface. By completely elucidating the bulk 
oscillations, this result paves the way for distinguishing the two-dimensional surface state in angular- 
dependent MR studies in Bi2Se3 with much lower carrier density. Besides, the present result provides 
a compelling demonstration of how the Landau quantization of an anisotropic three-dimensional 
Fermi surface can give rise to pronounced angular-dependent MR oscillations. 

PACS numbers: 71.18.-(-y, 73.25.-|-i, 72.20.-i, 72.80.Jc 



I. INTRODUCTION 

The three-dimensional (3D) topological insulator is a 
rapidly growing field of research in the condensed mat- 
ter physicSfi"— Because the bulk of a topological insula- 
tor belongs to a different Z2 topological classA^-^. than 
the vacuum, on the surface of a 3D topological insu- 
lator emerges an intrinsically metallic two-dimensional 
(2D) state which hosts helically spin-polarized Dirac 
fermionsiii^ Besides the profound implications of those 
spin-helical 2D Dirac fermions on future spintronics, the 
topological insulator is expected to present various ex- 
otic quantum phenomena associated with its non-trivial 
topology.—"— 

Motivated by theoretical predictions ^^lii three materi- 
als have so far been experimentally confirmedi^"— to be 
3D topological insulators: Bii_a;Sb2;, Bi2Se3, and Bi2Te3. 
The angle-resolved photoemission experiments l^^^iS have 
played decisive roles in indentifying the topologically 
nontrivial nature of their surface states. More recently, 
scanning tunneling spectroscopy experiment s ""^^ ' have 
elucidated the protection of the surface state from spin- 
nonconserving scattering. In addition to those surface- 
sensitive probes, transport experiments are obviously im- 
portant for understanding the macroscopic properties of 
topological insulators and for exploiting the applications 
of their novel surface state. However, capturing a sig- 
nature of the surface state of a topological insulator in 
its transport properties has proved difficult: The sur- 
face state has been seen by quantum oscillations only in 
Bii_a;Sbj; at a particular Sb concentration of 9%j2i the 
universal conductance ffuctuations associated with the 
surface state were observed only in Ca-doped Bi2Se3 af- 
ter careful tuning of the carrier densityi^ the Aharonov- 
Bohm oscillations through the surface state were ob- 
served only in very narrow nanoribbons of Bi2Se3 3^ Most 
often, the transport properties are dominated by the bulk 
conductivity due to residual carriers i^^i^ 

In this context, we have very recently found^S that the 
2D surface state of the topological insulator Bio.91Sbo.09 



gives rise to new types of oscillatory phenomena in the 
angular-dependence of the magnetoresistance (MR). It 
was demonstrated^^ that the oscillations observed at 
lower fields provide a way to distinguish the 2D Fermi 
surface (FS) in Bio.91Sbo.09, while another type of oscilla- 
tions of unknown origin was found to become prominent 
at higher fields. In the present work, we have tried to elu- 
cidate whether this new tool, the angular-dependence of 
the MR, can distinguish the 2D surface state in Bi2Se3, 
which is believed to be the most promising topological 
insulator for investigating the novel topological effects 
because of its simple surface-state structure j -^^'-*^^ We ob- 
served pronounced angular-dependent oscillations of the 
MR in an n-type Bi2Se3 single crystal with the carrier 
density rig = 5x10^^ cm"^, but our detailed analysis 
clarified that the observed oscillations are solely due to 
the Landau quantization of the anisotropic bulk FS of 
this material. Nevertheless, to the best of our knowl- 
edge, such an angular-dependent MR oscillations in a 3D 
material has never been explicitly demonstrated in the 
literature, so the present results nicely supplement our 
general understanding of the angular-dependent MR os- 
cillation phenomena. In addition, our analysis provides 
a solid ground for discriminating the contributions of the 
2D and 3D FSs in the angular-dependence of the MR in 
future studies of Bi2Se3 single crystals with much lower 
carrier density. 



II. EXPERIMENTAL DETAILS 

High-quality Bi2Se3 single crystals were grown by 
melting stoichiometric mixtures of 99.9999% purity Bi 
and 99.999% purity Se elements in sealed evacuated 
quartz tubes. After slow cooling from the melting point 
down to about 550°C over two days, crystals were kept at 
this temperature for several days and then were furnace- 
cooled to room temperature. The obtained crystals are 
easily cleaved and reveal a flat shiny surface. The X- 
ray diffraction measurements confirmed the rhombohe- 
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FIG. 1: (Color online) Temperature dependence of p^x in T. 
Upper inset shows the X-ray diffraction pattern of the Bi2Se3 
single crystal used for transport measurements. Lower inset 
shows Pyx for B \\ C3 measured at 1.5 K. The slope of pyx{B), 
shown by the thin solid line, suggests that the main carriers 
are electrons whose density is 3.4x10^* cm"''. 



dral crystal structure of Bi2Se3. Post-growth annealings 
at various temperatures under controlled selenium par- 
tial pressures were used to reduce the number of sele- 
nium vacancies that are responsible for creating electron 
carriers 

The resistivity p^^ was measured by a standard four- 
probe method on rectangular samples, with the electric 
current / directed along the Ci axis. Continuous rota- 
tions of the sample in constant magnetic field B were used 
to measure the angular dependence of the transverse MR 
within the C3-C2 plane. For some selected magnetic-field 
directions, the field dependences of pxx and the Hall re- 
sistivity Pyx were also measured by sweeping B between 
±14 T. 



III. RESULTS AND DISCUSSIONS 

A. Resistivity and SdH Oscillations 

Figure 1 shows the temperature dependence of pxx of 
the Bi2Se3 single crystal studied in this work. It shows 
a metallic behavior dp/dT > down to ~30 K, and 
saturate at lower temperature (there is actually a weak 
minimum near 30 K, as is usually observe d^^'^" in low- 
carrier-density Bi2Se3). The single-crystal nature of the 
sample is evident from the X-ray diffraction data shown 
in the upper inset of Fig. 1. The lower inset of Fig. 1 
shows the Hall resistivity pyx measured at 1.5 K for the 
field direction along the C3 axis, which suggests that the 
main carries are electrons and the carrier density rte is 
3.4x10^* cm~^ (in a one band model). From the values 
of the Hall coefficient Rh = 1-82 cm^/C and pxx = 0.28 
mJlcm at 1.5 K, the Hall mobility p,H is estimated to be 
6500 cmVVs. 
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FIG. 2: (Color online) SdH oscillations measured within the 
C3 — C2 plane. Thin solid lines are the result of our pxx{B) 
simulation (see text). Inset shows the measurement configu- 
ration. 



Figure 2 shows Pxx{B) measured at 1.5 K for several 
field directions in the transverse geometry (/ || Ci and 
B L I), after removing the antisymmetric components 
due to the leakage of pyx- Two features are evident: 
First, pronounced Shubnikov-de Haas (SdH) oscillations 
are seen for any field direction, suggesting their 3D origin. 
Second, the background of the SdH oscillations varies 
significantly with the field direction, indicating that the 
transverse MR is very anisotropic. Both features are 
taken into account in our simulation of the observed an- 
gular dependence of the MR, as will be discussed below. 

Figure 3 presents the analysis of the observed SdH os- 
cillations. The oscillations in dpxx/dB plotted as a func- 
tion of 1/B for B II C2 (0 = 90°) are shown in Fig. 3(a) 
as an example. The very simple pattern seen in Fig. 3(a) 
is a result of the single frequency F = 107 T (see inset for 
the Fourier transform) governing the SdH oscillations 1^ 
The same analysis was applied to the data for other field 
directions, and the obtained -F as a function of 9 are 
shown in Fig. 3(b). The same set of frequencies can be 
extracted from the Landau-level "fan diagram" [inset of 
Fig. 3(b)], which is a plot of the positions of maxima 
in Pxx{B) as a function of the Landau level numbers. 
The slopes of the straight lines in the fan diagram give 
exactly the same F{9) as the Fourier transform result. 
Another piece of information that can be extracted from 
the fan diagram is the phase of the oscillations, 7, which 
is determined by pxx ^ cos[2ti{F / B + 7)]. In the present 
case, all the straight lines in the inset of Fig. 3(b) in- 
tersect the horizontal axis at the same point, giving 7 = 
0.4 that is independent of the field direction. The an- 
gular dependence of the SdH frequency F{6) points to a 
single ellipsoidal FS located at the F point with the semi- 
axes ka = kb = 4.5x10^ cm-i (± C3) and = 7.3x10^ 
cm-i (II C3). The expected F{9) for this FS is shown by 
the solid line in Fig. 3(b), which fits the data very well. 
The carrier density corresponding to this FS is 5x10^^ 
cm~'^, which is about 50% higher than that obtained from 
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FIG. 3: (Color online) Analyses of the SdH oscillations. The 
inset in the middle schematically shows the obtained 3D Fermi 
surface and the definition of 0. (a) SdEI oscillations for B \\ C2 
as a function of 1/B. The Fourier transform shown in the in- 
set reveals a single frequency F — 107 T. (b) F{9) measured 
within the C3 — C2 plane; inset shows the fan diagram for sev- 
eral 6. (c) Temperature dependence of the SdH oscillations 
for B II C3; inset shows the temperature dependences of the 
SdH amplitudes measured along the C3 and C2 axes at 12 
T, yielding the cyclotron mass of 0.14me and 0.24me, respec- 
tively, (d) Dingle plots for B \\ C3 at several temperatures 
give the same To = 9.5 K. 



Rh, as was reported previously in the literature.— This 
discrepancy is due, at least partly, to the fact that the 
present Rh was measured in the low-field limit and con- 
tains the so-called Hall factor, which is usually between 
1 - 2. 

The temperature dependence of the SdH oscillations 
was measured for two field directions, along the C3 and 
C2 axes. Figure 3(c) shows the Pxx{B) data in B \\ C3 for 
some selected temperatures, where one can see that the 
background MR is essentially temperature-independent 
and that oscillations are still visible even at 25 K. The 
inset of Fig. 3 (c) shows the temperature dependence 
of the SdH amplitude measured at 12 T for B along C3 
(6* 0°) and C2 {9 ==90°). The fits with the standard 
Lifshitz-Kosevich theory'^'^ yield the cyclotron mass rric of 
0.14me and 0.24me for 9 = 0° and 90°, respectively. The 
energy dispersion near the conduction band minimum 
of Bi2Se3 is known to be parabolic, ■^^ and the observed 
anisotropy in TOc for the given FS is entirely consistent 
with the parabolic dispersion to within 5%. Relying on 
this dispersion, the Fermi energy Ep is calculated to be 
56 meV and the Fermi velocity vp is 3.8x10'' cm/s for 
any direction. 

The Dingle plots [shown in Fig. 3(d) for S || C3 as an 
example] yield the Dingle temperature Td = 9.5 K which 
is almost isotropic and is constant at low temperature; 
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FIG. 4; (Color online) (a) Angular dependences of pxx mea- 
sured within the C3 — C2 plane in constant magnetic fields, 
whose values were 0, 0.5, 1.5, 3, 4.5, 6, 8, 10, 12, and 14 T. 
(b) Simulation of pxx based on Eq. (1) for the same magnetic 
fields as in (a). 



this gives an isotropic scattering time t = 1.3x10""'^^ s, 
implying the mean free path £ (= vft) of ~50 nm and 
the mobility /isdH of about 1600 cm^/Vs, which is almost 
four times smaller than pLu estimated from Rh- Such a 
discrepancy between /isdH and pin has long been noted in 
various systemsj^^ and the essential reason lies in the dif- 
ference in how the effective scattering time comes about 
from the microscopic scattering process:"^^ The scattering 
rate l/rtr determining transport properties acquires the 
additional factor 1— cos0 upon spatial averaging (0 is the 
scattering angle), while 1/r to govern the dephasing in 
the quantum oscillations is given by a simple spatial av- 
eraging without such a factor. Hence, if the small-angle 
scattering becomes dominant (which is often the case at 
low temperature), l/rtr can be much smaller than l/r. 
This means that the mean free path in our sample can 
be even larger than that estimated above. 



B. Angular-Dependent MR Oscillations 

The new observation of the present work is the oscilla- 
tory angular-dependence of the MR shown in Fig. 4(a), 
where, on top of the two-fold-symmetric background an- 
gular dependence, pronounced oscillations are evident 
at higher fields. Importantly, the peak positions shift 
with magnetic field, which is different from the ordinary 
angular-dependent MR oscillations (AMRO) in quasi- 
low-dimensional systemsi^ In the following, we show 
that these oscillations are due to the Landau quantiza- 
tion of the 3D FS and, hence, are essentially of the same 
origin as the SdH oscillations. The SdH oscillations oc- 
cur as the Landau-quantized cylinders in the Brillouin 
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FIG. 5: (Color online) Fitting of the background MR for B \\ 
C2. Open circles show the experimental data [symmetrical 
component of p:cx{B)]; the dashed line demonstrates the low- 
field quadratic dependence of the resistivity, pqd(-B) = a2B^\ 
the dash-dotted line represents the high-field saturation limit, 
psat; the solid line is the resulting reconstructed background, 
Pbg{B) (see text). The left and right insets show the obtained 
angular dependences of the coefficients a2 and psat within the 
C3 — C2 plane. 



zone expand and cross the FS with increasing magnetic 
field, while the angular-dependent oscillations occur as 
the axis of those cylinders rotates in the Brillouin zone. 
Obviously, no oscillation is expected for rotating mag- 
netic field when the FS is spherical, but when the FS 
is anisotropic, the number of cylinders residing within 
the FS changes as the cylinder axis is rotated, leading to 
resistivity oscillations. Therefore, an anisotropy in the 
FS is a necessary ingredient for the angular-dependent 
oscillations. 

To show that the observed angular-dependent oscilla- 
tions are essentially due to the Landau quantization, it is 
most convincing to reconstruct the angular- dependence 
of the MR based on the SdH oscillations data. For this 
purpose, one needs to understand the exact magnetic- 
field dependence of the MR in the presence of the SdH 
effect, and know how it evolves when the magnetic field 
is rotated. Figure 5 shows the Pxx{B) data for B \\ C2 
[9 — 90°), which we take as an example for presenting 
our procedure to extract the necessary information. At 
low magnetic fields, the MR in Fig. 5 exhibits an al- 
most quadratic field dependence, which can be fitted by 
Pqd{B) = a2-B^. With increasing magnetic field, the MR 
tends to saturate, while pronounced SdH oscillations de- 
velop at the same time. The overall background MR 
Pbg{B) can be described as 1/pbg = 1/Pqd + l/Psat, 
which combines the low-field quadratic behavior and the 
high-field saturation of Pxx{B). This background devel- 
ops on top of the zero-field resistivity PxxiO) of about 
0.28 mJ7cm. The same fitting of the background of MR 
can be made for the whole range of the magnetic-field 
direction. The evolutions of the parameters 02 and psat 
with 9 in the C3 — C2 plane are shown in the left and right 
insets of Fig. 5. Note that the center of the symmetry 



for both 02 (9) and psat (9) is shifted from = 0° by about 
12°, which is probably due to the rhombohedral symme- 
try of the crystal which leads to the appearance of cross 
terms^^ in the magnetic-field expansion of the resistivity 
tensor 1^ 

Once pbg is known as functions of both 9 and B, one 
can superpose the SdH oscillations to obtain the expected 
Pxx{9,B). In the case of our Bi2Se3, thanks to the fact 
that the SdH oscillations are composed of only one fre- 
quency, one can simulate experimental data as 



Pxxi9,B) = pbg{0,B) X 

I + ARtRdRs cos 



27r 



m 

B 



+ 7 



(1) 



where Rt, Rd and Rs are temperature. Dingle, and spin 
damping factors, respectively 1^ As shown in Fig. 2 by 
thin solid lines, the experimentally observed SdH oscilla- 
tions are reproduced reasonably well with Eq. (I) with 
a parameter A ■ Rs = 0.12 which is independent of 9. 

The same set of parameters can be used to calculate 
the angular-dependent MR for fixed B, Pxx{9), with Eq. 
(1). The result is shown in Fig. 4(b), where one can 
easily see that the calculated angular dependences of the 
MR follow very closely the rather complicated patterns of 
the observed pxx (9) . This gives compelling evidence that 
the angular-dependent oscillations are essentially due to 
the Landau quantization of the anisotropic 3D FS. 

It was emphasized in Ref. [2^ that a merit of mea- 
suring the angular-dependence of the MR is that at suf- 
ficiently high magnetic field, the 3D FS remains in the 
quantum limit (where all the electrons condense into the 
1st Landau level) and do not contribute to the angular- 
dependent oscillations, while the 2D FS will always pro- 
duce MR oscillations when the magnetic-field direction 
is nearly parallel to the 2D plane. Once the density of 
the residual carriers in a topological insulator sample is 
sufficiently reduced, the 3D FS should be easily brought 
into the quantum limit and the 2D FS would become dis- 
tinguishable in the angular-dependent MR oscillations. 
Also, the SdH effect needs a certain range of magnetic 
field for the oscillations to be recognized, whereas the 
observation of the angular-dependent oscillations can be 
made only at the highest field; this gives a certain prac- 
tical advantage to the latter method when the quanti- 
zation condition is only barely achieved with the avail- 
able magnetic field. Therefore, one would expect that 
the angular-dependence of the MR will be a useful tool 
for distinguishing the 2D surface state of the topological 
insulator Bi2Se3 when a single crystal with much lower 
carrier density becomes available. 



IV. CONCLUSIONS 

In conclusion, we observed pronounced angular- 
dependent oscillations of the MR in high-quality single 
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crystals of n-type 612863 with rie = 5x10"'^* cm"^. We 
show, by simulating the angular-dependent oscillations 
based on the information obtained from the SdH analy- 
sis, that the oscillations are essentially due to the Landau 
quantization of the 3D bulk Fermi surface. This provides 
a compelling demonstration of how the Landau quantiza- 
tion of an anisotropic 3D FS can give rise to pronounced 
angular-dependent MR oscillations. Furthermore, the 
present results pave the way for distinguishing the 2D 
surface state in 612863 in future studies of the angular- 



dependent MR, by completely elucidating the oscillations 
due to the bulk FS. 
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